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ABSTRACT 



Aims. We present a framework to obtain photometric redshifts (photo-zs) for gamma-ray burst afterglows. Using multi-band photom- 
etry from GROND and SwiftfUVCfT, photo-zs are derived for five GRBs for which spectroscopic redshifts are not available. 
Methods. We use UV/optical/NIR data and synthetic photometry based on afterglow observations and theory to derive the photo- 
metric redshifts of GRBs and their accuracy. Taking into account the afterglow synchrotron emission properties, we investigate the 
application of photometry to derive redshifts in a theoretical range between z ~ 1 to z ~ 12. 

Results. The measurement of photo-zs for GRB afterglows provides a quick, robust and reliable determination of the distance scale 
to the burst, particularly in those cases where spectroscopic observations in the optical/NIR range cannot be obtained. Given a suf- 
ficiently bright and mildly reddened afterglow, the relative photo-z accuracy 77 = Az/(1 + z) is better than 10% between z = 1.5 and 
z ~ 7 and better than 5% between z = 2 and z = 6. We detail the approach on 5 sources without spectroscopic redshifts observed with 
UVOT on-board Swift and/or GROND. The distance scale to those same afterglows is measured to be z = 4.3LjJ"[j for GRB 080825B, 



2.13!" " for GRB 080906, z 



3.44+^? for GRB 081228, z: 



2.03^ ;° for GRB 081230 and z 



1.28^ !° for GRB 090530. 



Conclusions. Due to the exceptional luminosity and simple continuum spectrum of GRB afterglows, photometric redshifts can be 
obtained to an accuracy as good as 77 ~ 0.03 over a large redshift range including robust (77 ~ 0.1) measurements in the ultra-high 
redshift regime (z > 7). Combining the response from UVOT with ground-based observatories and in particular GROND operating in 
the optical/NIR wavelength regime, reliable photo-zs can be obtained from z ~ 1.0 out to z ~ 10, and possibly even at higher redshifts 
in some favorable cases, provided that these GRBs exist, are localized quickly, have sufficiently bright afterglows and are not heavily 
obscured. 

Key words, gamma rays: bursts, Techniques: photometric 
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1. Introduction 

The measurement of cosmological redshifts using photome- 
try in broad-band filters avails of prominent and character- 
istic features in the spectra of extragalactic objects, such as 
the strong 4 000 A Ca H/K break in early-type galaxies (e.g. 
B au m1ll962[). or the Lvm an-bre aks for highly redshifted objects 
(e.g.lSteidel & Hamiltonlll992l:ISteidel et al.lll996r,lMadau et alj 
1996). The crude spectral coverage of broad-band photome- 
try compared to even low-resolution spectroscopy limits both 
the accuracy and the application of photo-z measurements to a 
redshift range where the data contain significant redshift sig- 
natures. However, with the advent of large-field, deep multi- 
band photomet ric and spectroscopic surveys as the SPSS (e.g. 
YorketalJl2000h . the Hubble Deep Fields ( e.g. IWilliams etal] 
19961). COSMOS fe.g . l^viUeet alJl200l . or GOODS (e.g. 



Giavalisco et alj2004l) two major limitations of photometric red- 



length coverage and a lack of training sets of spectroscopic 
redshifts to calibrate and validate the photometric redshift de- 
terminations convincingly. Robust and accurate photo-zs can 
now be obtained for a large number of objects efficiently and 
to a depth which is ina ccessible to spectroscopy at even the 
largest telescopes (e.g. Sawicki e t al.l I l997t iHogg et alJ [1998; 
Csab aiet al l l2003t iFranxetalJ 120031: iMobasher et alj 2004; 
Caoak et alJl2007t lllbert et al J 120061 l2009t ISalvato et all [2009; 
Bouwens et al]l2010UMcLure et alJl201QY 



A further application for photometric redshifts in addi 
tion to multi-object and faint-source distance determination 
is i n the field of gamma-ray burst (GRB) astronomy (see 



I2009at lOates et al J [20091 



2009: 



e.g. Tagliaferri et a02005bljakobsson et al.ll2006bt Haislip et al 
2006; ICurran et al l 120081 : iRossi et al l l2008bt iGreiner et al 



Tanvir et al 

120091: lOlivares et al.1 12009at iMcBreen et all 1201 a Cucchiara 
et al., in prep.). Despite its extreme early luminosity, a 
GRB afterglow fades quickly. This enables an unambiguous 



Salv aterra et al 
20101 



shifts were overcome: limited photometric accuracy and wave- 
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identification but requires that afterglow spectroscopy (e.g . 



opy 

Jakobssonetal] l2006bb iKawai et al.l T2006: 



Berger et al.l 120051: 
Prochaska et al . 2007a; Fvnbo et al.ll2009l) must be performed in 
a timely manner in the first few hours to days after the prompt 
y-ray emission. Rapid multi-band observations of the afterglow 
with small to medium aperture telescopes provide photometric 
redshift measurements in cases where spectroscopic observa- 
tions do not reveal emission/absorption lines, or are unfeasible, 
for example when the optical/NIR afterglow is either too faint, 
has not been identified to arcsec accuracy or the spectroscopic 
resource is unavailable. Photometric follow-up on timescales 
of several tens of seconds to minutes and the resulting photo-z 
also provides the trigger for a finetuned setup for l ong-slit spec- 
trographs with hmited wavelength coverage (e.g. iGreiner et alj 
2009b: Ra tiltani20Toh . 

The weal th of photometric r edshif t codes (e.g . 

Connolly et all 1 1995b fBolzonella et alJ l2000t iBenftezI l20Qi 




00; 
lav 



2-10' 
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Bender et al. 2001; Tagliafe rri et al.l 120021: iCollister & Laha . 
2004 Carlil es et a can essentially be sub-divided into 



Wavelength [A] 
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learning based methods and template fitting. The earlier requires 
a sufficiently large training set of photometry together with 
spectroscopic redshifts, which although limits the application 
(object types, different redshift/magnitude space), is however 
self-contained and makes no prior assumptions on the physical 
properties of the unknown objects. In contrast, template methods 
assume a spectral shape as obtained from observations and 
models of similar objects and compare its synthetic photometry 
with observations of the source of interest, which extrapolates 
reasonably well into unmeasured redshift ranges. The probabil- 
ity of a redshift solution is then evaluated using, for example,^ 2 
or Bayesian statistics. 

The natural approach for GRB afterglows is template based, 
as it incorporates the well-known emission properties of the light 
source. The emission of GRB afterglows arises when the ultra- 
relativistic ejecta from the GRB central e ngine are decelerated 
by the swept-up circumburst medium (e.g.[P aczvnski & Rhoads 
19931; Meszar os & Reesl 119971: [PirarJ 120041: iMeszarosI 12006: 
20071: IGehrels et alj 120091) 



Fig. 1. UVOT (uvw2, uvm2, wvwl, u, b, v) effective areas 
(against left y-axis) and GROND (g', r 1 , /', z', J, H, K s ) fil- 
ter curves (right y-axis), respectively in colored areas. The 
UVOT effective areas have been obtained from the most recent 
HEASARC Calibration Database 1 . GROND filter curves include 
all optical components in GROND including the telescope, but 
exclude the atmosphere. Shown in black lines are template af- 
terglow spectra, from redshifts z=l, 2, 3, 4.5, 6 to 8 (Top left 
to bottom right). These spectra also differ in their spectral index 
and rest-frame extinction (amount and reddening law), see e.g. 
the two synthetic spectra at z=4.5 and 6, which show a redshifted 
2 1 75 A dust feature. Also the DLAs centered at 1 2 1 .5 x ( 1 +z) nm 
differ in their hydrogen column densities. 

v~^t~ a , and concordance (Q« = 0.27, Oa = 0.73, Hq 
7 1 km/s/Mpc) cosmology. All errors are given at lcr confidence 
unless indicated otherwise. 



leased via synchrotron emission from external shocks, where 
the kinetic particl e energy is transformed into radiation (e.g. 
Wiie rs et alj|1997l) . Hence, the theoretical continuum spectrum 
of afterglows is a three-fold broken power law, with breaks 
at the self-absorption frequency v fl , the inj ection or t ypical 
frequency v,- and the cooling frequency v c ( Sari et al.l [1998; 



The optical afterglow is re- 2. Photometric Redshifts for GRB afterglows 



iPanaitescu & Kumarll2000l: iGranot & Sarill2002h . In the late af- 
terglow, v, < v c (the slow cooling regime) and the optical 
wavelength range is located either above or below v c , which is 
fully supported by current observations (e.g. iGalama & Wiiers] 
2001: lStratta et alJ2004l:lNardini et alJ2006r.lSchadv et al.l2007t 
Greiner et alj 120101) . Two prominent signatures, the Ly-a and 



Lyman-limit break at 121.5 x (1 + z) nm and 91.2 x (1 + z) nm 
respectively, and their characteristic redshift-dependent promi- 
nence, allow for a robust and precise redshift determination. 

Here we present a framework of GRB afterglow photo-z 
measurements, investigate their uncertainties, and detail the ap- 
proach on several afterglows observed with the Swift/Ultra Violet 
Optic al Telescope (UVOT; IGehrels et all l2004t Ikoming et all 
120051) covering the 170-600 nm wavelength range in six fil- 
ters and th e Gamma-Ray burs t Optic al/Nearlnfrared Detector 
(GROND; IGreiner et all 120071, 120081) which is sensitive be- 
tween 360 nm and 2300 nm using seven filter bands 
(g', r>, i', z', J, H, K s ). 

Throughout the paper we adopt the convention that the flux 
density of the afterglow F v (v, t) can be described as F v (v, t) oc 



2. 1 . Filters and Photometry 

The key ingredient for any photometric redshift measurement 
is high-quality photometry over a large wavelength range. For 
GRB afterglows Swift/UWOT and GROND offer the natural 
data source as both instruments systematically follow-up on 
GRB triggers and nicely complement each others sensitivity 
and wavelength coverage. UVOT onboard the Swift satellite is a 
30 cm space-based telescope primarily sensitive in the ultravio- 
let (UV) and optical range using uvw2, uvml, uvwl, u, b, v fil- 
ters, that starts o bserving the GRB fi eld as quickly as ~40 s 
after the trigger (Romi ng et al.l I2009T) . Additional optical and 
near-infrared (NIR) response is provided by GROND, a seven 
channel imager (g', r', i', z', J, H, K s simultaneously) mounted 
at the 2.2 m MPI/ESO telescope at the ESO/LaSilla observatory. 
GROND's reaction to GRB triggers is hence subject to visibil- 
ity constraints, and is typically in the range of a couple to sev- 
eral hours post-burst. The effective area of the UVOT filtersC] and 
GROND's sensitivity curves are all shown in Fig.Q] 



2.2. General constraints on afterglow photo-zs 

The photometric identification of a redshift signature in the spec- 
tral energy distribution (SED) of a GRB afterglow requires de- 



http://heasarc.gsfc.nasa.gov/docs/heasarc/caldb/ 
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tections in at least two filters to measure the continuum and an- 
other one to locate the Lyman-break. In this minimal case, the 
continuum is not well constrained, of course, and a strong degen- 
eracy exists between afterglow spectral index and intrinsic red- 
dening. Consequently, the redshift information is rather crude. 
Additional observations in independent wavelength ranges pro- 
vide an improved continuum determination or constraints on the 
spectral index, and result in a more accurate photo-z determi- 
nation. The condition of having at least three individual filters 
might possibly be relaxed if further information about the spec- 
tral index or the dust content of the afterglow is availabl e from 
X-ray measurements as provided by Swift/XKT dBurrows et al 
2005) and the optical to X-ray flux ratio (e.g. lAfonso et al 
2010). In the standard model, for example, the optical spectral 



dening similar as observed in the SMC (e.g. Stratta et all 2004; 
iKann et al.ll2006LlStarling et alJ2007tlSchadv et al.ll2007ll2010l) . 
several recent observations show a strong 2175 A bump in the 
afterg low's SED at the redshift of the burst (e.g. iKruhler et al 
2008; IProchaska et"ai] 120091: lElfasdottir et al.ll2009h iDeng et al 



index {S is required to be p o - fix if both wavelength ranges 
probe the same part of the synchrotron spectrum , or (3 n _ 0X-O.5 
in case of a cooling break lying in between (e.g. Sari et al. 1998; 
iGranot & Sari 2002]). Also, a large Ay OSt solutions could be con- 
sidered unlikely in cases where the combined fit requires the de- 
reddened optical to X-ray flux ratio (3 x > fix- This is when the 
extrapolation of the dust-corrected UV/optical/NIR SED would 
significantly over-predict the X-ray measurement. Similar com- 
bined optical/X-ray SED analysis make however strong assump- 
tions on the nature of the emission in X-ray and optical energy 
ranges, in particular that both are emitted by the same popu- 
lation of radiating electrons. The differ ences in the respective 
light curves (e.g. Panaitescu et al.ll2006l) . and t heir inconsisten- 
cies with the most simpl e fireball scenarios (e.g.lWillingale 'etafl 
120071; [Evans et al.1 12009) raises questions, whether this is indeed 
the case for all afterglows. Therefore, the most reliable approach 
is solely based on UV/optical/NIR data. 

Even in the case of a well-sampled SED, several systematic 
uncertainties limit the accuracy of the photometric redshift de- 
termination, where the most dominant ones are the uncertainty 
in the dust-reddening properties, in the opacity of the Ly-a for- 
est and the strength of the Damped Lyman-o' Absorber (DLA) 
associated with the GRB. To quantify the systematic effects 
and the accuracy of the redshift measurement via photometry, 
a large sample of 4000 GRB afterglow SEDs was simulated as 
described in Sec. l2.3l and analyzed following Sec. 12.41 



2.3. A synthetic afterglow spectral energy distribution 

The continuum emission of an afterglow in the UV/optical/NIR 
range is synchrotron radiation, and usually described as a sin- 
gle power law {F V (A) = F(){A/Ao~f, for caveats on this as- 
sumption, see Sec. 12. 9\ . This continuum spectrum is affected 
by reddening due to dust in the host galaxy. Therefore, dif- 
ferent dust columns with average attenuation laws as observed 
in the Milky Way (MW), Small and Large Magellanic Clouds 
(SMC, LMC) in the parametrization according to Hi3 d 19921) 
between the burst and the observer are added to the synthetic 
spectrum. In addition, a more generic approach is used to ob- 
tain a broader range of plausible dust exti nction laws via th e 
Drude model proposed in lLi et alJ (12008b and lLiang &H d2009). 
Local dust-reddening laws strongly differ in their absolute UV 
absorption, which for a given Ay decreases in strength from 
SMC, over LMC to MW. The most prominent extinction fea- 
ture is the 2175 A bump, which is generally attributed to ab- 
sorpti on by graphite grains (e.g. S tec her & Pond [19651 : iDraind 
120031) . While the feature is highly significant in MW and LMC 
models, it is absent in the SMC dust attenuation law. Although 
most bright afterglows are best described with featureless red- 



2009: iPerlev et alJ l2010b). where its location can be used as a 
redshift tracer in case the photometric data are well sampled. 
A large number of optical afterglow spectra also show 

the 



presence o f a stron g DLA (e .g. Jensen et ail 2001 
;swii k et alJ 12004 IChenet alJ 120051: IWatsonetal 



Vreeswiik et al. 2004; Chen etalJ 120051: I Wa tson et al. 
20061: ISavaglia 120061: iFvnbo et alJ 120061: iBerger et alTfe oOf: , 
Jakobsson etail l2006al: IProchaska e~ al. 2007bt lFvnboetal.1 
2009). To quantify the effects of neutral hydrogen absorption 
associated with the GRB host, DLAs with different hydrogen 
columns are added to the after glow spectrum following the 
description of iTotani et alJ (l2006h and references therein. 

Bluewards of the DLA centered at the redshifted A a , i.e. 
at 121.5(1 + z) nm in the observers fram e, the after glow flux 
is further suppressed by the Ly-a forest dLvndsl IT97lh : inter- 
vening neutral hydrogen absorbers between the burst site and 
the observer. To account for different measurements, measure- 
ment errors and specific sight lines, the Ly-a effective opti- 
cal depth TLy- ff is al lowed to vary according to th e constraints 
given in Table 4 of iFaucher-Giguere et al.l d2008l) in the red- 
shift range between 2 and 4.2. Below z ~ 2, TL y -„ is essen- 
tially zero, and above z ~ 4.2, T\^- a and its error are extrap- 
olated from the lower-redshift data. The derived optical depth 
TLyc is then converted into an averaged flu x depression factor 
(D) dOke & Korvcanskvlll982HMadaulll995l) . 

Combining the upper effects, the synthetic afterglow spec- 
trum becomes (see Fig.[T]for examples): 



F v (A,x) = (D a (z))F ^— j exp[-T dust (z,A v )-T DLA (z,N H )] (1) 

for Ap(l + z) < /l bs < &a(J + z) and similar for {Dp) and higher 
order hydrogen absorptions. Here, x are the free parameters 
x = {z, Ay,P,N H )- Below the Lyman-limit at ~ 91.2(1 + z) nm, 
we assume that the observed flux is fully attenuated by the neu- 
tral hydrogen along the line of sight, and hence F V (A) = 0. The 
intrinsic brightness term Fq of the afterglows is selected to be in 
a range of previously observed UVOT/GROND afterglows and 
accounts for the typical reaction time of the instruments. 

Synthetic AB magnitudes of the afterglow in the different 
filters i are derived via: 

f A- l F v (A,x)T(A)dA 

mag -=- 2 - 51og j^vm +219mag (2) 

where T(A) are the specific filter curves as shown in Fig. [T] 
and F v is given in yuJy. The synthetic afterglow magnitudes are 
then varied using a Gaussian probability distribution with the 
associated photometric errors as its standard deviation. These 
errors are the superposition of a constant, i.e. minimum error 
according to the absolute photometric accuracy obtained with 
standard calibration and a brightness dependent term related to 
photon statistics. The earlier is conservatively set to 0.04 mag 
for mvw2, uv m.2, uvwl, u, b, v and,?', r', i', z' an d 0.05 mag for 
/, H, K s (cp. lPoole et alJ2008tlGreiner et alJ[2 008). The latter is 
a function of the brightness of the simulated afterglow, and takes 
into account the sensitivities of UVOT and GROND as well as 
their reaction time to GRB triggers. In detail, the brightness of 
the simulated afterglow has been mapped to a magnitude error 
using the image statistics of standard follow-up observations (2 h 
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Fig. 2. Relation between brightness of an object and photomet- 
ric accuracy detailed for a 2 h GROND and 600 s UVOT ex- 
posure under fairly typical observing conditions (3 days from 
new moon, airmass 1.5, 1 arcsec seeing). The black dots show 
the brightness and photon noise of the individual objects in the 
r'-band, while the solid line denotes the relation between ob- 
ject brightness and total (photon+calibration) magnitude error. 
Dotted (for GROND) and dashed (for UVOT) lines show the 
same relation for the remaining filters. Individual objects are 
not plotted for these filters to enhance clarity. Horizontal dashed 
lines indicate typical limiting magnitudes for 2cr, 3<x and 5<x 
confidence levels. 



GROND integration at 8 h after the trigger, 10 min UVOT inte- 
gration at 1 h after the trigger) as shown in Fig. [2] 

The resulting magnitudes and associated errors are then used 
to derive a photo-z as described in Section l2~4l to compare in- and 
output values. 

2.4. photo-z code 

The software used to derive afterglow photometric redshifts i s 
based on the publicly available hyperZ (Bo lzonella et al.l l2000). 
which minimizes the^ 2 from synthetic photometry of a template 
spectrum F V {A, x) against the observed data, i.e.: 



F y.obsU;) - const, x F V (A, x) 



(3) 



where Fy^iAf) denotes the measurements in different broad 
band filters i with associated errors cr, . 

In addition to the exis ti ng tre atment of the Lyman absorp- 
tion a ccording to iMadaul d 1995b and default reddening tern- 



plates (Allen 1976; Seaton 1979; Fitzpatrick 1986; Prevotet al 
1984; B ouchet et alj 119851; ICardelli et al.1 119891; ICalzetti etal 
2000), the code is complemented by several additions. In par- 



ticular, power-law spectral templates with the possibility to con- 
strain the spectral ind ex, a DLA and a reddening law following 
Maio lino et alJ d2004l) were added. It also includes the measured 
total UVOT and GROND filter response, including all optical 
components in the pathway from the primary mirror of the tele- 
scope JoJh£_c|Ujmtun^fficj^ (see Fig. Q] and 
also lGreiner et aUfeOOSfclPoole et alj|2008l) . 

The results are similar to the standard outputs of the original 
hyperZ version, including the photometric redshift, error bars or 
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Fig. 3. Properties of the sample of 4000 simulated GRB after- 
glow spectra with respect to their host extinction Ay OSt , spectral 
index f3, simulated redshift z s ; m , and DLA hydrogen column den- 
sity log(A^n) from top left to bottom right. 



contours of the three main parameters (redshift, spectral index 
and intrinsic reddening) at arbitrary confidence levels and the 
statistical probability associated with the derived redshifts and 
secondary solutions when relevant. 

2.5. Properties of the afterglow mock sample 

The properties of the mock set of 4000 afterglow spectra are 
shown in Fig. [3] The sample properties are chosen to be as close 
as possible to what is known abo ut optical aftergl ows with re- 
spect to their spectral indices (e.g. Kan n et alJ2010h and the neu- 
tral h ydrogen column densities of their DLAs (e.g. iFvnbo et"ail 
2009) but are not fully representative of the global GRB after- 
glow properties. All previous demographic studies of GRB after- 
glows are still strongly biased against highly-reddened and red- 
shifted bursts, and complete afterglow properties are hence sub- 
ject to large uncertainties. Instead, the simulated sample is used 
for the determination of typical accuracies and systematic effects 
in the photo-z measurement, when the afterglow is detected by 
UVOT and/or GRO ND. In part i cular, the dust distribution w hich 
is based on iKann et alj d2010h , ISchadv et ail d2007l 1201 Oh and 
iGreiner et al.l d2010h is probably heavily skewed towards low 
dust environments. It hence rather resembles the sight lines to- 
wards bright and mildly extinguis hed afterglows as typically de- 
tected in the UV/optical r ange dKann et alJl2006t ISchadv et al.l 
2007; Ro ming et al]l2009l) . The redshift distribution is chosen to 
peak at the redshift interval where mo st of the Swift bursts origi- 
nate from (z ~ 1 -4. lFvnbo et alJ20 09). but includes a significant 
number of ultra-high redshift GRBs (~1000 at z > 8) to inves- 
tigate the application of photometric redshifts at these extreme 
values. 



2.6. The application and accuracy of afterglow photo-zs 

Figure [4] shows a comparison between redshifts of previous af- 
terglows obtained from spectroscopy against photo-zs derived 
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from GROND and UVOT data (see also Tab. Q]) and summa- 
rizes the result of the simulation. In the simulation, all available 
extinction laws (see Sec. 12.4-b are used in the analysis, and the 
reported photo-z is based on the reddening law that returns the 
best-fit SED, i.e. the minimum^- 2 . 

Below z ~ 0.8, even UV photometric measurements do not 
cover the wavelength range of the Lyman-limit absorption, and 
the redshift is hence unconstrained. The sensitivity towards low- 
redshift events begins at z ~ 1 .0, when the Lyman-limit cuts off 
a significant amount of flux in the uvw2 filter. The relative uncer- 
tainty rj = Az/(1 +z) in the low -redshift regime z < 1.5, however, 
is still comparatively large with r\ > 0.08. In addition, there is a 
significant fraction of afterglows (~ 20-40%), where the photo- 
metric data do not provide redshift constraints (Az < -0.5). The 
observations do not cover the part of the spectrum blue-wards 
of the spectral break, and Ly-limit absorption is somewhat de- 
generate to intrinsic dust extinction for UV detections with low 
S/N and large associated photometric errors. The redshift sig- 
nature of the Ly-break can hence not be reliably identified in 
these cases. This effect is clearly visible in Fig. [4] as an under- 
prediction of the average redshift of the sample by photometric 
measurements in this redshift interval. 

The dust-redshift degeneracy is broken with increasing red- 
shift, as the Ly-limit moves to redder wavelengths, producing a 
drop-out in uvw2 at z > 1.7 and in uvm2 at z > 2.0, which is 
too sharp to be mimicked by dust. Accordingly, 77 decreases to 
77 ~ 0.05 over this redshift range. At higher redshifts, the ab- 
solute redshift uncertainty stays roughly constant and hence the 
relative scatter decreases to 77 ~ 0.04 at z ~ 4 until redshifts of 
z ~ 5, when it further drops to 77 < 0.03. At these redshifts, ab- 
sorption in the Lyman-a forest becomes significant (see Fig. [T}, 
and hence there are two signatures in the spectrum which are 
used to measure the redshift. Furthermore, there are enough fil- 
ters both blue-wards, and red-wards of the breaks to measure 
their location and the continuum with high accuracy. 

At even higher redshift, the accuracy remains essentially 
constant until z ~ 6.5. This demonstrates that the total numbers 
of individual filters does not strongly affect the robustness of the 
photo-z measurement, as long as the continuum is fairly well de- 
termined. The number of filter bands that contain constraining 
information decreases from 13 at z ~ 2 to 5 at z ~ 6.5. Due to 
the intrinsic power-law spectrum of GRB afterglows, just a few 
filters red-wards of the break are sufficient to reliably measure 
the continuum, at least in the case of blue (/3 < 0.5, Ay < 0.2 
with 2 filters) or mildly red (J3 < 1.2, A v < 0.4 with 3 filters) 
events. Information about the location of either Ly-o- or the Ly- 
limit is then derived via 1 or 2 filters covering the wavelength 
range at or blue-wards of the break(s). 

As the wavelength spacing between z' and J is relatively 
large, the photometric redshift is rather loosely constrained be- 
tween z ~ 6.5 and z ~ 8.3, increasing from 77 ~ 0.04 or Az ~ 0.3 
at z = 6.5 to 77 ~ 0.1 or Az ~ 1 at a redshift of around 8. Also at 
redshifts z > 8, photometric colors can put strong constraints on 
the redshift. If the source is relatively bright as compared to the 
sensitivity limit (as input in the simulations), the H - K s color, 
coupled with the flux decrease in J and the z'-band upper limit, 
yields a robust redshift measurement (77 ~ 0.1) until z ~ 9.5. 
There are, however, a significant number of rather strong outliers 
with Az > 1.5 in this region, which are associated with moderate 
dust reddening along the line of sight. Above a redshift of z ~ 10, 
the transmitted flux in J is below the instruments sensitivity limit 
even for intrinsically bright afterglows. With detections in only 
H and K s , the measurement is fully degenerate between the spec- 
tral power-law slope /3, and in particular redshift and reddening. 



Only if the /-band limit is deep enough to exclude dust as the 
cause of the extremely red SED, a rough (around Az ~ 1-2) 
redshift estimate might be obtained under specific assumptions 
and possibly in combination with X-ray observations of the fad- 
ing afterglow until redshifts of around z ~ 12. 

The accuracy of photometric redshifts is further demon- 
strated and verified by photo-zs obtained for a number of after- 
glows where spectroscopic redshifts could be secured (red data 
in Fig. H]i. These sources span a redshift range from z - 1.67 to 
z ~ 8.2, and support the results of the simulation. 

2.7. Systematic effects 

In the redshift range above z > 6.5, two systematic effects bias 
the measurement (blue lines in Fig. 0). In contrast to the opti- 
cal/UV filters, the NIR filters do not overlap due to emissivity of 
the NIR sky in certain wavelength ranges. In cases where the red- 
shift signature is located in the filter gaps, it can not be measured 
precisely of course. For example, photometric measurements 
yield very similar SEDs for z ~ 7 and z ~ 8 afterglows with 
only upper (z < 8.3) and lower (z > 6.5) limits on the photomet- 
ric redshift, and there is hence a trend of over/under-estimating 
the redshift around z ~ 6.5 - 7 and z ~ 7.5 - 8, respectively. 
A T-band filter centered at 1020 nm, as used in the an alysis 
for GRB 090423 dTanvir et alJ2009tlSalvaterr"aet al. 2009), pro- 
vides additional information in this redshift range. Furthermore, 
the continuum is only sparsely-probed at these redshifts. With 
three detections in JHK S , different combinations of little dust, 
blue /3 , and a lower-z as well as no dust, red /3 and a higher red- 
shift fit the input data equally well. While the returned best-fit 
photo-z is based on the solution with lowest possible dust con- 
tent, the input Ay OSt values follow Fig. 12.31 Hence, the average 
photo-z of the sample consequently somewhat overpredicts the 
simulation input value, but with errors containing the input pa- 
rameter space, of course. The redshift dependence of the average 
photo-z (thick line in Fig. |4j> is the net-result of both aforemen- 
tioned effects. 

The asymmetry of the blue-shaded error regions in the two 
lower panels of Fig. [4] can be readily understood as the re- 
sult of the upper systematic effects. The best-fit average photo- 
z is based on the no-dust solution, but a family of solutions 
with somewhat higher dust values, which as a consequence also 
means lower redshifts, describe the simulated measurements 
equally well. This is represented by the asymmetric error regions 
extending to lower redshifts in the lower panels of Fig.|4]and also 
illustrated in Fig. [5] 

Furthermore, secondary and possibly tertiary solutions at 
lower redshifts and larger Ay OSt s of around 1-5 mag exist for 
afterglows where the continuum is constrained by only two fil- 
ters, which includes all ultra high-z events (z > 8.5). Given a 
sufficiently bright afterglow, and hence accurate photometry, the 
low-redshift solutions can generally be ruled out at >99% con- 
fidence under the assumption of a conventional dust attenuation 
law. 

2.8. The effect of photometric accuracy 

The absolute accuracy in the photometric measurement is the 
basic quantity in the accuracy of photometric redshifts (see also 
e.g. iBolzonella et alJl2000t lllbert et al.l 12009b . Also the break- 
ing of the dust-redshift degeneracy strongly depends on pho- 
tometric accuracy in cases where the continuum spectrum is 
constrained by only few filters. Fig. [5] shows the z - A^ ost con- 
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Fig. 4. Photometric redshift accuracy. Small black dots show the 
mock set of simulated afterglow spectra and their corresponding 
photo-z. Thick blue lines show the average photometric redshift 
after distributing the 4000 mock afterglows into redshift bins of 
100 afterglows each. The lowest two panels also show in blue- 
shaded areas the quadratic sum of the typical difference to the in- 
put redshift and the lcr statistical uncertainty of the photo-z anal- 
ysis averaged over 100 afterglows in absolute (Az = z p h t - Zsim) 
as well as relative (rj = Az/(1 + z)) terms. The statistical prob- 
ability that the true value is in the interval of +/-20% around 
the quoted lcr accuracy is >95%. Grey shaded areas represent 
the Zphot > constraint. The thin, horizontal, dotted lines in 
the lowest two panels represent Az =-1, -0.5, 0, +0.5, +1 and 
77 = Az/(1 + z) =-0.1, -0.05 -=0, +0.05 and +0.10. The large 
red dots show final UVOT/GROND photo-z measurements for 
real bursts where a spectroscopic redshift has been obtained (see 
Tab.Q]). The green dot shows the photo-z of the flat-spectrum ra- 
dio quasar PKS 0537-286 derived in a s imilar manner (z = 3.10; 
iBottacini et aIll2O10Hwri"ght et al.lll978l) . 



Table 1. Photo-zs of GRB afterglows compared against spec- 
troscopic redshifts . Errors on spectroscopic measurements are 
reported when relevant (> 0.01) and available. 



GRB z (spectroscopic) 



070802 
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080804 
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090516 
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lutions. The redshift is no longer constrained, and the 3<x contour 
nearly contains the full redshift interval. 



tours of GRB 09 0423 (z ~ 8.2 - 8.3. ISalvaterraeTall l2009l 
ITanvir et"aDl20 09) for different fixed photometric errors. Here, 
we used g', r' , /', z', J, H, K s measurements only (excluding 
the available T-band imaging) for a direct comparison of the re- 
sults. Two aspects are clearly apparent. Firstly, and already men- 
tioned in Sec. 12.61 and shown in Fig. [4] is the asymmetric shape 
of the contours including regions of lower redshift with some 
dust contribution. Secondly, the large increase in the allowed pa- 
rameter space when going from highly accurate (+0.03 mag) to 
crude (±0.25 mag) photometric measurements. In the latter case, 
photometry can no longer disentangle high-z and large A^ ost so- 



2.9. Caveats in the photo-z measurement 

The key uncertainty in the photometric redshift measurements 
is the attenuation of photons due to dust along the line of sight, 
which is generally assumed to be similar to local galaxies or 
to previously observed extinction laws. The dust extinction at 
high redshift, and in particular in extreme environments such 
as the circumburst medium of a GRB might be significantly 
different. The dust around the GRB might possibly be subject 
to destruction by the intense UV radiation of t h e afterglow or 
the GRB progenitor dWaxman & Draind l2000t iDraine & Had 
I2002t iD'Elia etalJ 120071) . albeit a clear and highly significant 
observational signature of this process is still lacking (e.g. 



T. Rriihler et al.: Afterglow photo-zs 



7 




Fig. 5. Photometric redshift accuracy for GRB 090423 with dif- 
ferent photometric errors. From top left to bottom right, the 
analysis wa s perfo rmed with the g', r', i', z', J, H, K s data of 
iTanvir et a l. (2009) but with photometric errors of +0.03 mag, 
+0.06 mag, +0.15 mag and +0.25 mag respectively. The black 
line corresponds to the l<x contour, while the increasingly dark 
shaded areas correspond to the 90%, 95.4% (2cr), 99%, 99.73% 
(3cr) and 99.99% confidence contours. The vertical dashed line 
denotes the spectroscopic redshift. 



iPerlev et al.l201 0a). Although nearly all afterglo w SEDs are well 
described with continuous extinction laws (e.g. ISavaglio & Fall 



. . _ .. s (c , 

2004 IStratta et alJ)2004t iKann et al]l200fil IStarling et alJl2007t 



Schadv et alj I2007L 1201 Ol) . a very sharp break in the attenua- 
tion curve could possibly be misidentified as the signature of the 
Lyman breaks. This dust feature, however, must ev en be sharper 
than what was observed for GRB 070318 dWatsonll200"9l) , so far 
the only case where such a break has been reported. 

A second caveat is present in the case of sources located 
at regions with high G alactic foreground re ddenings. The avail- 
able foreground maps ( Schle geTet alJ[l9 98') are limited to a spa- 
tial resolution of few arcminutes, and can hence be subject of 
uncertainties in the foreground correction of up to several tens 
of percent. For GRBs, typically at large Galactic latitudes with 
Eb-v 5= 0.1 this effect is usually smaller than the absolute photo- 
metric accuracies and hence negligible. Extreme caution is how- 
ever required when interpreting spectral breaks of objects close 
to the Galactic plane (| b |< 5 - 10°) associated with regions of 
Eb-v Z 0.5. A detailed analysis is crucial in these cases. 

Finally, there is the possibility that a break in the syn- 
chrotron spectrum, most likely the cooling break at v c is lo- 
cated in or evolves through the UV/optical/NIR range at the 
time of the observations. To quantify the effect of an evolv- 
ing cooling frequency through the optical bands, the previously 
used continuum emission of a single power-la w was replaced by 
a smoothly connected broke n power-law (e.g. lBeuermann et al.l 
Il999t iGranot & Sari] 120021) . The difference between the two 
power-law slopes /3\ - @2 and the smooth ness of the break were 
set to 0.5, and 1, respectively, following Grano t & Saril (2002) 
for a cooling break in the slow-cooling regime. Fig. [6] shows 
photometric redshifts obtained for four mock afterglows with 
a fixed flux at 1000 nm restframe of 60 /j.Jy, where the break 
wavelength evolves through the UV to the NIR. These synthetic 
afterglows are located at redshifts of z = 1.5, z = 2.5, z = 3.5, 
and z - 4.5, and hence span the range of the largest fraction 
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Fig. 6. Photometric redshift accuracy for four synthetic bursts 
(z = 1.5, z — 2.5, z — 3.5 and z - 4.5) with a smoothly broken 
power law as continuum spectrum against the location of the 
break wavelength. Horizontal dashed lines mark the input red- 
shift, and the best-fit photo-z for an unbroken continuum, and 
black solid lines the average best fit photometric redshift when 
following Sec. 12.41 Grey shaded areas represent the typical l<x 
uncertainty, derived in similar manner as in Sec. 12.61 



of Swift GRBs starting from the low redshift end of the photo-z 
sensitivity. Their photo-z with respect to break wavelength and 
compared to the standard single power law continuum spectrum 
is shown in Fig. [6] The introduced curvature due to the cooling 
break is interpreted as an increased dust content, but can not re- 
semble the strong breaks due to Lyman-blanketing. Hence, even 
in case of a broken power-law continuum spectrum, the photo- 
metric redshift can be considered reliable. 

2. 1 0. Number statistics 

The total number of bursts where a photo-z could in princi- 
ple be derived with the presented method depends on a num- 
ber of factors, primarily a UVOT or GROND detection and 
wavelength coverage of the redshift signatu re. Using the typical 
GROND and UVO T detection efficiency (iRoming 'et all [2009; 
iGreiner et al.ll20i p j) and t he most recent GRB redshift distri- 
bution dFvnbo et al.ll2.009h . and requesting a favorable declina- 
tion for GROND observations and a Galactic latitude cut of 
| b |> 10°, the presented method is applicable to around 15-30% 
of all GRBs detected by SwiftfBAT, and 30-50% of all GRBs 
with a detected optical afterglow. For many (70%-90%) of the 
latter sources spectroscopic observations will of course be fea- 
sible and return a much more accurate redshift measurement. 
Nevertheless, in few of all bursts (~ 1-5%), which tend to include 
the rare LAT or ultra high-z events (e.g. IGreiner et alJ l2009a. 
Cucchiara et al., in prep.) a photo-z provides the best redshift 
measurement for GRBs and their afterglows. 

3. Application to individual bursts 

To further demonstrate the concept of photometric redshifts for 
GRBs, a number of afterglows between July 2007, which is 
when GROND began systematic follow-up observations of all 
GRBs, until May 2010 were extracted from the UVOT and 
GROND archives. These afterglows must be located at regions 
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of low Galactic foreground reddenings {Eg-v < 0.3), detected 
by UVOT/GROND in at least three filters, and no spectroscopic 
redshift for the bursts is available. Finally, the SED of the after- 
glows must show a prominent break in the observed wavelength 
range to constrain the redshift to at least 90% confidence. Five 
out of a total of ~ 105 GROND observed bursts (i.e. around 5%) 
meet the earlier constraints. Out of these five sources, four trig- 
gered Swift/BAT and one triggered SuperAGILE. 

3.1. Data reduction and cross calibration 

GROND data for the selected bursts shown in Tab. [2] were re- 
duced in the standard ma nner using pyraf/IRAF dTodvl 119931: 
iKiipcu-Yoldas et all [2008). The data in different filters are ob- 
tained simultaneously by hardware setup with the exception of 
a small difference in the mean photon arrival time betwee n the 
g', r', /', z' and JHK S filters (see e.g. lKruhler et alJl2009bl) . At 
the average observing time of GROND after the GRB trigger of 
several hours, this difference is 5f/f b. s < 10~ 3 , and hence negligi- 
ble in the analysis. Absolute photometry fo rg', r', i', z' has b een 
tied to the SDSS standard star network (Smit h et alj|2002l) or 
nearb y fields covered by the SDSS catalog DR7 (Abazaiian et alj 
2009). Photometry fo r JHK S has been deri ved against 2MASS 
field stars in all cases dSkrutskie et alj|2006l) . The GROND mea- 
surements, and the respective calibration source are given in 
Tab. E] 

Swift/UVOT photometry has been obtained following 
iPoole et al] d2008l) . As UVOT operates through filter cycles, the 
measurem ents had to be int erpolated to a common epoch de- 
tailed in Scha dv et al.l d2010t) . Typically, the interpolation is per- 
formed over a short time interval, and the introduced uncertain- 
ties are much smaller than the individual measurement errors 
dominated by photon noise. UVOT photometry and the refer- 
ence time of its observations are provided in Tab. [3] 

The reference time of UVOT measurements is generally 
much earlier than the ground-based observations, but the large 
overlap in the bv and g'r' filters (see Fig. Q]) offers a straight- 
forward cross calibration. Using the synthetic photometry from 
the mock spectra in Sec. [2] and including stellar templates, color 
terms are derived over a large range of photometric colors: 

b-g' =0.15(g'-r') + 0.03(g'-r') 2 V(g'-r')e [-1,2] (4) 

V-r' = 0.62(r'-/') +0.10(r'-/') 2 V(r' -/') e [-1,2] (5) 

g' - b = -0.20(6 - v) - 0.05(6 - vf V (b - v) e [- 1 , 2] (6) 

r' - v = -0.55(6 - v) - 0.04(6 - vf V (6 - v) e [- 1, 2] (7) 

where all UVOT magnitudes are in the AB system. The system- 
atic error on the color terms is of order 2% for the individual 
equations, and smaller than 5% for all of them. Due to the large 
spectral overlap of the GROND and UVOT filters, the system- 
atic uncertainties introduced through the cross calibration pro- 
cess are between 5 and 10 mmag even for SEDs with red colors 
of g' — r' ~ 1. This is much smaller than the individual mea- 
surement errors (see Tab. and does not introduce additional 
uncertainty in the measurement. 

If possible, the transformation equation with the largest 
spectral overlap (i.e., Eq. |4|i is used to transform the UVOT 
measurements to the GROND epoch, where errors are prop- 
agated accordingly. This implicitly assumes that the over- 
all spectral evolution of the optical transient associated with 
the burst is achromatic. While there is evidence from well- 
sampled optical afterglow light-curves of chromatic evolution 



Table 4. Properties of the GRB afterglows in the photo-z sample 
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simultaneous to strong variability (e.g. iGreiner et al.l l2009d) 
the associated color changes are generally mode rate to absent 
even i n complex panchrom a tic light-curves (e.g . Masetti et al. 



2000; Holland et al. 2003; Covino et al. 
2004[ kacusin et al.l 120081: lupdike et al l 
2008[|Bloom et al.|l2009l:IOates et al.l2009l:lKruhler et al 
Cenko et al.l l2009£ iPerlev et all l2010at iMcBreen et al 



Linki n et al 



Perlev et al 



2009a; 



2010; 



iKann et al . 2010; C ovino et al.l2010l) . Hence, even if a chromatic 
light-curve evolution is present, either due to the passage of the 
cooling break, or associated with late inner engine activity, it is 
very unlikely to mimic the strong breaks in the UV/optical/NIR 
SED due to absorption by neutral hydrogen used for the photo- 
Z measurement. In all cases, however, the hydrogen column of 
the DLA is unconstrained by the photometric measurement, and 
a DLA with log Nh ~ 21.5, which is about the median value 
for afterglows wh ere such a measurement was possible (e.g. 
Fvnb o et alj2009t) is adopted in the SEDs (Figures0[8]|9][IO]and 
fTTb . The properties of the five afterglows including their photo- 
metric redshifts are summarized in Tab. [4] 



3.2. GRB 080825B 



GRB 080825B wa s detected by AGILE dEyangelista et aU 2008) 
and Konus-Wind dGolenetskii et al.l l2008|). and the optica l/X- 
ray afte r glow was rapidly identified by iThone et all d2008l) and 
iPaganil §008). The g\ r' , i', z', J, H, K s GROND SED ob- 
tained at a midtime of 7.2 h after the trigger (Fig. |7]i shows the 
clear presence of two strong breaks between g' and r' and r 1 
and i', respectively (see Fig. [7]). In addition there is evidence 
for curvature red-wards of 650 nm which is well described with 
a SMC-type reddening. The best-fit Of 2 = 1.68 for 3 d.o.f) 
photometric redshift is z p h t 

of A^, ost = 0.20 + 0.15 and a spectral index of $ 
consi stent with the redshift limits from Ide Ugarte Postigo etafl 
(200 8J). At this redshif t, and using the y-ray properties from 
Golenetskii et al. (2008), the isotropic equivalent energy release 
of the GRB in the rest-frame 1 keV to 10 MeV energy range is 
x 3.7 x 10 53 erg with a rest-frame peak energy of ~ 380 keV. 

3.3. GRB 080906 



4.31+qIj with a host extinction 
A4+0.3 



GRB 080906 dVetere et al.ll2008l) trigg ered Swift, and its optica l 
afterglow was detected by UVOT 82 s dHolland & Veterdl2008l) 
and with GROND 10.7 h dAfonso et al.ll2008bl) after the burst. 
The combined SED extends from u to the K s band and is shown 
in Fig. [8] The photometric redshift is driven by the uvwl non- 
detection and is z p h t = 2.1 3^20 with negligible intrinsic host 
extin ction and a spe ctral index of ft = 0.5+y 2- Previous claims 
from lHollandl f2008) are consistent with this redshift. 
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[magvej 


[magvesa] 




080825B 


7.20 


22.45 ± 0.07 


19.91 : 


tO.04 


18.62 ±0.04 


18.24 ±0.04 


16.83 ± 0.04 


16.03 ±0.05 


15.30 ±0.05 


SA103-626/2M 


080906 


10.9 


22.11 ±0.07 


21.72: 


t0.05 


21.49 ±0.06 


21.37 ±0.07 


20.18 ±0.16 


19.55 ± 0.23 


18.75 ±0.32 


SDSS/2M 


081228 


0.719 


22.10 + 0.14 


20.79 : 


tO.04 


20.30 ± 0.05 


20.12 ±0.08 


18.64 ±0.14 


17.63 ±0.14 


16.87 ±0.15 


SDSS/2M 


081230 


5.44 


21.62 ±0.04 


21.16: 


t0.03 


20.98 ± 0.04 


19.75 ±0.04 


19.62 ± 0.05 


18.99 ±0.06 


18.50 ±0.14 


SA94-242/2M 


090530 


21.8 


22.15 ±0.14 


22.08 : 


h0.10 


21.78 ±0.14 


21.64 ±0.18 


> 20.33 


> 19.50 


> 18.74 


SDSS/2M 



(fl) All magnitudes are observed magnitudes, i.e. uncorrected for the corresponding Galactic foreground reddening. 

m Calibration sou rce for the data. All magnitudes from SA standard stars and fields are taken from the primary S loan standard star network 
dSmith et alJl2002h . 2M and SDSS d enote a calibration a gainst stars from the 2MASS catalog dSkrutskie et all 120061) in case of J, H, and K s 
measurements, and the SDSS DR7 ( Abazajian et al. 2009) for g', r', i', z' band data. 



Table 3. UVOT photometric data 



GRB 


Af <fl) 


uvw2 (b) 


uvm2 


uvwl 


u 


b 


V 




[h] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


080906 
081230 
090530 


0.194 
0.139 
1.67 


> 19.54 

> 19.82 

21 17+0.62 
ila -0.39 


> 19.01 

> 19.62 

19 74+0 25 

ly - '^-0.20 


> 20.16 

> 19.76 

20.29 +0 ) « 


1 Q Q5+0.40 

19 31+ 030 


> 18.90 

1926 +o.i9 

20 l q+8:18 


18 93 +u 11 

19 61+ 073 

20 36 +0 48 



UVOT magnitudes are the result of interpolation to the given reference time. 
All data are observed magnitudes, and in the UVOT system. 




No reddening 



3, 




Fig. 7. Broad-band spectral energy distribution of the afterglow 
of GRB 080825B as observed with GROND. The foreground- 
corrected SED shows two prominent breaks corresponding to 
Lyman-limit absorption in the g' band, and Lyman-o- in r' . The 
corresponding redshift is z p h t = 4.31+jj [ji. The inset shows the 
confidence contours of the redshift solution versus intrinsic ex- 
tinction, where the three-dimensional z - /3 - A^ ost parameter 
space has been collapsed onto a two-dimensional z - A!j, ost grid. 
The increasingly dark shaded areas correspond to the 68.3% 
(lcr), 90%, 95.4% {2(f), 99%, 99.73% {3(f) and 99.99% con- 
fidence contours, where the 68.3%, 90% and 99% contours are 
also marked with solid lines. The different filter bands are plot- 
ted with their effective wavelength, and the horizontal error bars 
represent their FWHM, i.e. from 50% to 50% of maximum trans- 
mission. 



3.4. GRB 081228 



Swift/BAT triggered on GRB 081228 (lPageetal.1 12 008), 
and GROND detected the optical afterglow (Afons o ;t alj 
2008a), while UVOT observations only yield upper limits 



Fig. 8. Foreground corrected broad-band spectral energy distri- 
bution of the afterglow of GRB 080906 with UVOT (open cir- 
cles) and GROND (filled circles). Upper limits are shown by 
downward triangles. The SED shows a strong break due to the 
Lyman-limit being located in the uvwl band. The photometric 
redshift is z p h t = 2. 1 3^^'2q- The inset shows the confidence con- 
tours of the redshift solution versus spectral index, where the 
three-dimensional z - - Ay OSt parameter space has been col- 
lapsed onto a two-dimensional z — /? grid. Lines and shadings are 
the same as in Fig. [7] 



(lLandsman & Pagell2008l) . The SED in the GR OND filters is 
domin ated by a red g' - r' color of ~ 1.3 mag dAfonso et al.l 
2008c), which yields a photometric redshift of z p h t = 3.44^5 "ij. 
There is no strong evidence of excess absorption with a best fit 
spectral index f3 



0.03+j} ( '3 assuming 



1.24+01' with A host 
a SMC type reddening law (x 2 = 2.34 for 3 d.o.f). The SED is 
slightly better fit with a 2175A feature (Zphot = 3.49+° ^ and x 1 = 
1.58 for MW-like reddening for 3 d.o.f, and z phot = 345+° ^ and 
X 1 =2.14 for LMC-like reddening for 3 d.o.f, see Fig. [9j, which, 
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Fig. 9. Foreground corrected broad-band spectral energy dis- 
tribution of the afterglow of GRB 081228 as observed with 
GROND. The SED shows a prominent break corresponding 
to Lyman-a in the g' band. The resulting redshift is z p hot = 



o 44+0.15 

-0.32- 



Fig. 10. Foreground corrected broad-band spectral energy dis- 
tribution of the afterglow of GRB 081230 as observed with 
UVOT (open circles) and GROND (filled circles). Upper limits 
are shown as downward triangles The SED shows a prominent 
break corresponding to Lyman-limit absorption between uvwl 
and u bands. The resulting redshift is z p hot = 2.03^ [4. 



given the small improvement inx 2 , is only significant at the 0.9 
and 0.5cr level, respectively. 

3.5. GRB 081230 

The UV/optica l /NIR afterglow of the Swift GRB 081230 
(LaParolaetal. 2008) was imaged by UVOT starting 131 s 
(bates & La Parolall2008l) and GROND 4.2 h after the trigger 
(Ufonso et alJ \2QQ% . The combined SED (Fig. [TOj is broad, 
ranging from the u to K s filter, whereas the afterglow is unde- 
tected in the UV bands. The implied break between the uvwl 
and u bands yields a photometric redshift of z p hot = 2.03*2 }^. 
There is mild curvature in the SED which is well fit (x 1 = 8.2 
for 9 d.o.f) with a moderate amount of SMC-type reddening with 



an A h v ost = 0.22 



+°-J* and a spectral index /3 = OA^ 2 ,. LMC and 



MW-like reddening models are strongly ruled out due to the ab- 
sence of a 2175A dust feature and are not shown in Fig.fTOl 

3.6. GRB 090530 



GRB 090530 was detected by Swift (ICannizzo et alJ 120091) . 
and its afterglow was observed with space- and a number 
of ground-based teles copes (e. g. iNissinen & Hentunenl 1 2009: 
iFlewelling etail 12009]) . UVO T riSchadv & Cannizzoll2009h~ and 
GROND dRossietal.ll2009cT) observations started 2.5 min and 
21.3 h after the bursts. The combined UVOT/GROND SED is 
shown in Fig. [TT] and contains all filters except the NIR, which 
given the late GROND observations only yield non-constraining 
upper limits. The data are acceptably fit (y 2 = 1 A for 6 d.o.f.) 
with a marginally {A* ost = 0.15+ 015 
law of spectral index /3 = 0.4 
models provide slightly worse fits to the data, but within the er- 
rors comparable redshifts, extinctions, and spectral indices. The 
dominating spectral feature is a break between the two bluest 
UVOT filters, which implies a photometric redshift of z p h ot = 



og ) SMC-type reddened power 
: 0.3. LMC and MW extinction 



1 to+0.15 „ 

-0.17' Z ph ot 



i To+0.17 
10 -0.16' Z P hot 



1.37+0 is for SMC ' LMC and 



MW extinction laws, respectively. Lower-redshift (z < 1) solu- 
tions are allowed at the 2<x level. 



SMC reddening 
LMC rcddeiiiiii 
MW reddening 




Fig. 11. Foreground corrected broad-band spectral energy distri- 
bution of the afterglow of GRB 090530 as observed with UVOT 
(open circles) and GROND (filled circles). The SED exhibits a 
break indicative of Lyman-limit absorption between the uvwl 



and uvml bands. The corresponding redshift is z p h t = 1.28^ [ 6 



4. Conclusions 

GRB afterglow photometry in multiple bands offers a viable and 
robust method to derive the distance scale to the burst to reason- 
able accuracy. In cases where the photometric observations cover 
the wavelength range blue- and redwards of the redshift tracers 
such as the Lyman-limit in case of low z < 5, or Lyman-ff for 
high (z > 3) redshift GRBs, photo-zs can be obtained in princi- 
ple for redshifts from z ~ 1 out to z ~ 12 for bright events with 
state-of-the-art follow up (e.g. Oates et al. 20091; [Bloom et al. 
120091: iTanvir et alJ [20091 ISalvaterra et al.l 12009V The accuracy 
of the photometric redshift determination is a function of red- 
shift, wavelength coverage and quality of the photometric mea- 
surement, and can reach relative errors as good as rj ~ 2 - 3% 
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for observations as routinely performed by systematic follow-up 
with SwiftfUVOT and GROND. 

The photo-z method presented in this work is based on stan- 
dard template fitting, and derives the best-fit solution of redshift, 
reddening and spectral index of the afterglow against the ob- 
tained data simultaneously. In this way, redshifts for five GRBs 
between z = 1.2 and z = 4.3 with relative errors rj < 0.1 and 
without spectroscopic observations are derived using photomet- 
ric UV/optical/NIR data from UVOT and GROND. 

Due to the simple synchrotron continuum, and the absence of 
stron g reddening observed for a good fracti on of all afterglows 
(e.g. iKann et aO l2010fc iGreiner et al.l | 2010|) and in particular 
the high-redshift ones (eTgJKann et al.ll2007l;lRuiz-Velasco et al ^ 



2007; Greineretal. 2009b; Tanviretal. 2009; Salvaterra et al 



2009: IZafaretal] l2010). two filters can constrain the spectrum 



such that a strong break due to Ly-a can reliably be identified in 
a third one. There is, however, a strong observational bias against 
highly-redshifted, highly-reddened afterglows, and a first mod- 
erately extinguished, z ~ 5 afterglow has only been detected re- 
cently dPerlev et al.ll2010ah . Follow-up observations in the stan- 
dard NIR broad-band filters can yield robust photo-z measure- 
ments for bright events up to distance scales where Lyman-ar is 
redshifted into the H band. Due to the uncertainty in the dust ex- 
tinction, around 10% or 30% of the well-detected high-z events 
beyond z = 9 are expected to be deviating from the true redshift 
with more than 15% or 10%, respectively. It must be stressed, 
again, that these results are obtained from simulated afterglows 
with pre-defined reddening properties and brightnesses in a way 
that they are amenable to detailed photometric studies. 

The strong breaks, their characteristic prominence and the 
well-defined continuum emission also provide a unique redshift 
solution under the assumption of a conventional dust attenuation 
law, and hence prevent catastrophic errors due to a degeneracy 
in spectral features such as a confusion between the Lyman-limit 
and the 4000 A break for galaxies fe.g. lSawicki et alJll997l) . 

Albeit with much larger uncertainties than spectroscopic red- 
shifts, the accuracy of afterglow photo-zs is adequate for fur- 
ther studi es about the aftergl ow reddening along the line of 
sight (e.g. .[Greiner et alJl2009ah . the soft X-ray ab sorption (e.g., 
Rossi et al.ll2008bl) , the GRBs energy budg et (e.g. lGreiner et alj 
2009ar. lAbdo et al]l2009t [Amati et alJ.120091). its emission mech- 
anisms (e.g. lAbdo et alJ 120091 : iKumar & Barniol Duranl l2009t 
IZhang & Pe'erll2009l) and the role of GRBs as probes for quan- 
tum g ravity ("e.g. lAmelino-Camelia & Smohnll2009t Ellis et al 



2009), the extragalactic background li ght (e.g. Gilmore et al. 
2009) or demog raphic studies of GRBs ( Jak obsson et al . 2006b; 
Fvnbo et al.l2009l) . Furthermore, photometric observations reach 



deeper flux limits than spectroscopy in particular in the NIR, 
which is of primary importance for redshifts z > 7. In case of 
these ultra-high redshift events, an identification through pho- 
tometry coupled with follow-up of their hosts with the upcom- 
ing generation of telescopes and instruments could reveal their 
hosts, and hence possibly a galaxy in the process of its first star 
formation. 
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